use "dataset_east_baton_rouge_cross", clear

forvalues year = 2000(10)2020 {

	cor votingblackpct registeredblackpct if year == `year'
	local myrho = string(r(rho), "%20.03fc")
	local mytitle = "`year' (r=`myrho')"

	#delimit;

	gr tw
		(sc votingblackpct registeredblackpct if year == `year',
			col(black) msym(Oh))
		(function y=x, range(0 1) lcol(black))
		,
			legend(off)
			xtitle("")
			ytitle("")
			subtitle(, fcol(none) lwid(none))
			plotregion(style(none))
			ylab(, angle(horiz))
			title(`mytitle')
			name(g`year', replace)
			;
			
	#delimit cr
	
	gr close

}

gr combine g2000 g2010 g2020, col(1) ///
	b1title("Share Black of those registered") ///
	l1title("Share Black of those voting") ysize(2) xsize(1)
	
gr export "_output/figure A6.pdf", replace

